All Fiber Vector Magnetometer Based on Nitrogen-Vacancy Center

Magnetometers based on nitrogen-vacancy (NV) centers in diamonds have promising applications in fields of living systems biology, condensed matter physics, and industry. This paper proposes a portable and flexible all-fiber NV center vector magnetometer by using fibers to substitute all conventional spatial optical elements, realizing laser excitation and fluorescence collection of micro-diamond with multi-mode fibers simultaneously and efficiently. An optical model is established to investigate multi-mode fiber interrogation of micro-diamond to estimate the optical performance of NV center system. A new analysis method is proposed to extract the magnitude and direction of the magnetic field, combining the morphology of the micro-diamond, thus realizing μm-scale vector magnetic field detection at the tip of the fiber probe. Experimental testing shows our fabricated magnetometer has a sensitivity of 0.73 nT/Hz1/2, demonstrating its feasibility and performance in comparison with conventional confocal NV center magnetometers. This research presents a robust and compact magnetic endoscopy and remote-magnetic measurement approach, which will substantially promote the practical application of magnetometers based on NV centers.


Introduction
The nitrogen-vacancy (NV) center consists of a substitutional nitrogen atom (N) combined with a vacancy (V) in one of the nearest neighbor sites of the diamond crystal lattice. Due to the excellent properties of long-lived spin coherence and fluorescence stability, nitrogen-vacancy centers in diamond offer an unprecedented solid-state platform for magnetometers and thermometers with hypersensitivity and ultra-high spatial resolution as an atomic-scale defect in the diamond crystal [1][2][3][4][5][6]. NV center magnetometers have the potential of miniaturization and practicality at room temperature under Earth's magnetic field background. Other state-of-the-art magnetometers, such as superconducting quantum interference devices and spin-exchange relaxation-free magnetometers, exhibit sensitivities at the fT/Hz 1/2 level, which is a thousand times higher than NV center magnetometer's sensitivity. However, they require more rigorous working conditions and longer measurement times, which greatly limit their application in practical engineering. In contrast, NV center magnetometers have the characteristics of fast start-up and can work in complex physical field environments such as high temperature and high pressure. NV center magnetometers have been widely used in room-temperature measurement of single living cells [7], high resolution imaging of neural activity [8], as well as magnetic field imaging of advanced material [9,10] in recent years.
However, the most excellent results of NV center magnetometers were achieved in a lab environment with a bulk device. It is necessary to develop compact, integrated, and mobile NV sensors for industrial applications. Fiber integration is an essential method to enable NV center sensors for practical application. Fiber-based NV center sensors have been providing compact and practical solid-state tools for bio-imaging, chip inspection, and nondestructive testing [11][12][13][14][15]. Some studies have been conducted by attaching diamonds varying in size from tens to hundreds of microns at the end of the multi-mode fiber, taper fiber, and photonic crystal fiber for highly efficient laser excitation and fluorescence collection, with magnetometer sensitivity in the ranges of 1.5-300 nT/Hz 1/2 [16][17][18][19][20]. Fiber taper and photonic crystal fiber were applied to enhance fluorescence collection efficiency because of their large numerical aperture. These fiber-based NV center sensors have also been expanded to magnetometer scanning probe or remote measuring probe. However, there still remain several space light optical elements that need to be precisely spaced, which makes the system unstable and inflexible. Moreover, the space light path is easily interfered with in practical use. Space optical elements, such as objective and dichroic lenses, greatly limit the practical application of the NV center magnetometer. In previous studies, a 1 × 2 fiber optical coupler was applied instead of a dichroic mirror to separate laser and fluorescence [20,21], but its relatively low optical transmission efficiency restricted the sensor's performance.
To address the aforementioned issues, we present an all-fiber NV center system based on multi-fiber optical elements with high fluorescence efficiency and extraordinary optical transmission efficiency. In addition, the all-fiber NV center magnetometer can be quickly set up outside the lab and operate in complex environments, which presented a robust and compact magnetic endoscopy and remote magnetic measurement approach. Due to the small size of the diamond particles, the all-fiber NV center magnetometer could also realize µm-scale magnetic field detection. Optical interrogation of the micro-diamond in the confocal system and fiber system was investigated to estimate the optical performance of the NV center system. The all-fiber NV center magnetometer was achieved with a sensitivity of 0.73 nT/Hz 1/2 , presenting a robust and compact magnetic endoscopy and remote magnetic measurement approach.
In some published works by other groups, a three-axes Helmholtz coil and repetitive optically detected magnetic resonance (ODMR) scans with magnetic field rotation were used to identify the relative orientation of NV center axes in the lab-frame coordinate and measure the direction and amplitude of the vector magnetic field [22,23]. This method was complicated and time-consuming. Based on the all-fiber NV magnetometers, we also proposed a simplified and rapid approach for vector magnetic field measurement, considering the relationship between the morphology of micro-diamond and NV centers orientations [24]. Using a static magnetic field with a known direction B a provided by one axis coil, we could determine the relative orientation of the NV center axes and establish the transformation relationship between lab-based Cartesian system and NV center axis coordinate system with a single ODMR scan. Following that, ODMR scans of the measured magnetic field B m and superimposed magnetic field B a +B m are conducted to measure the amplitude and direction of B m in lab-based Cartesian system.

Materials and Methods
The fiber-based magnetometer is shown in Figure 1a. The 532 nm laser (100 mW) is first sent to a fiber-connected acousto-optic modulator (AOM, S-M200-0.4C2C-3-F2S, Gooch & Housego, Ilminster Somerset, UK) to achieve a quick-switch laser pulse for pulsed ODMR measurement. Compared with continuous ODMR measurement, pulsed ODMR technique avoids optical and MW power broadening of the spin resonances, which can achieve higher sensitivity. Then the laser is guided to the wideband multi-mode circulator (WMC, WMC3L1F, Thorlabs, NJ, USA) from Port 1. About 98% of the laser power from Port 1 is transmitted to pump the diamond, which is glued at Port 2, as shown in Figure 1b. About 92% of the fluorescence collected by Port 2 (multi-mode fiber, a core = 200 µm, NA 1 = 0.39) largely spreads to Port 3. The light transmitted from Port 1 to 3 and Port 2 to 1 is isolated. The laser and fluorescence emitted from the diamond are effectively separated by the WMC, replacing the traditional dichroic mirror. The green and red arrows in Figure 1a show the path and direction of the laser and fluorescence, respectively. At the same time, about ten times higher optical transmission efficiency was achieved compared to 90:10, 1 × 2 fiber optical coupler. The fiber collimator, long-wavelength filter, and focusing lens were fixed to an avalanche photodetector (APD, 430A/M, Thorlabs, NJ, USA) through lens tubes and SM1 screw thread to collimate, filter, and focus fluorescence from Port 3 onto the detector's active area.
Nanomaterials 2023, 13, x FOR PEER REVIEW 3 of 12 laser power from Port 1 is transmitted to pump the diamond, which is glued at Port 2, as shown in Figure 1b. About 92% of the fluorescence collected by Port 2 (multi-mode fiber, a core = 200 μm, NA1 = 0.39) largely spreads to Port 3. The light transmitted from Port 1 to 3 and Port 2 to 1 is isolated. The laser and fluorescence emitted from the diamond are effectively separated by the WMC, replacing the traditional dichroic mirror. The green and red arrows in Figure 1a show the path and direction of the laser and fluorescence, respectively. At the same time, about ten times higher optical transmission efficiency was achieved compared to 90:10, 1 × 2 fiber optical coupler. The fiber collimator, long-wavelength filter, and focusing lens were fixed to an avalanche photodetector (APD, 430A/M, Thorlabs, NJ, USA) through lens tubes and SM1 screw thread to collimate, filter, and focus fluorescence from Port 3 onto the detector's active area. To manipulate the electron spin state of NV centers, a radio-frequency signal was generated by a signal generator (N5171B, Keysight, CA, USA). It was then routed through a microwave switch (ZFSWA2R-63DR-S+, Mini-circuits, NY, USA), an amplifier (ZHL-16W43-S+, Mini-circuits, NY, USA), and finally delivered to the copper coil antenna (Figure 1b). The oscilloscope performed the data acquisition of the detected signal from the APD. Both the AOM and the microwave switch were both controlled by a Programmable TTL Pulse Generator (PG, PulseBlasterESR-500 MHz, Spincore, FL, USA), which is mounted on the computer's motherboard. To identify the relative orientation of the NV center axes and measure the vector magnetic field, a one-axis Helmholtz coil was used to generate a biased magnetic field.
The micro-diamond we employed was a type-Ib high-pressure, high-temperature (HPHT) diamond with a diameter of about 550 μm, as shown in Figure 1c. The diamond was irradiated with an energy of 10 MeV and a dose of 1 × 10 18 cm −2 electron beam, then annealed at 850 °C in vacuum and 500 °C in the air, yielding a density of NV centers approximately 1.08 × 10 18 cm −3 . The micro-diamond was a cubo-octahedron with surfaces of 8 hexagons and 6 quadrilaterals in morphology. The hexagon surfaces correspond to the <111> crystal plane of the diamond. It is effortless to ensure thaqt the <111> crystal orientation of the diamond coincides with the axial direction of the optical fiber by gluing a hexagon surface to the fiber end face. The micro-diamond was deposited on the tip of Port 2 under an optical microscope by ethyl cyanoacrylate glue. To improve the mechanical strength of the fiber probe, the fiber end near the micro diamond was encapsulated by a ceramic ferrule. The diamond was exposed out of the ceramic ferrule, and the copper coil antenna was wound around the ceramic ferrule near the diamond. To manipulate the electron spin state of NV centers, a radio-frequency signal was generated by a signal generator (N5171B, Keysight, CA, USA). It was then routed through a microwave switch (ZFSWA2R-63DR-S+, Mini-circuits, NY, USA), an amplifier (ZHL-16W43-S+, Mini-circuits, NY, USA), and finally delivered to the copper coil antenna (Figure 1b). The oscilloscope performed the data acquisition of the detected signal from the APD. Both the AOM and the microwave switch were both controlled by a Programmable TTL Pulse Generator (PG, PulseBlasterESR-500 MHz, Spincore, FL, USA), which is mounted on the computer's motherboard. To identify the relative orientation of the NV center axes and measure the vector magnetic field, a one-axis Helmholtz coil was used to generate a biased magnetic field.
The micro-diamond we employed was a type-Ib high-pressure, high-temperature (HPHT) diamond with a diameter of about 550 µm, as shown in Figure 1c. The diamond was irradiated with an energy of 10 MeV and a dose of 1 × 10 18 cm −2 electron beam, then annealed at 850 • C in vacuum and 500 • C in the air, yielding a density of NV centers approximately 1.08 × 10 18 cm −3 . The micro-diamond was a cubo-octahedron with surfaces of 8 hexagons and 6 quadrilaterals in morphology. The hexagon surfaces correspond to the <111> crystal plane of the diamond. It is effortless to ensure thaqt the <111> crystal orientation of the diamond coincides with the axial direction of the optical fiber by gluing a hexagon surface to the fiber end face. The micro-diamond was deposited on the tip of Port 2 under an optical microscope by ethyl cyanoacrylate glue. To improve the mechanical strength of the fiber probe, the fiber end near the micro diamond was encapsulated by a ceramic ferrule. The diamond was exposed out of the ceramic ferrule, and the copper coil antenna was wound around the ceramic ferrule near the diamond.

Optical Interrogation of NV Centers in Diamond
The all-fiber NV center magnetometer was analyzed in terms of optical properties and sensitivity performance, compared with the traditional confocal system. Laser loss Nanomaterials 2023, 13, 949 4 of 11 in the fiber system mainly occurred at AOM, fiber mating sleeves(FMS), and WMC. The diffraction efficiency of fiber-coupled AOM was 70%. Multi-mode fiber connectors provided low insertion loss, typically less than 0.5 dB, and the insertion loss from single-mode to multi-mode was negligible. The laser power transmitted to the diamond reached 64.2 mW when the laser source power was set to 100 mW, as measured by an optical power meter. In addition, the laser transmission efficiency of the traditional space light confocal system was about 62.3%, as shown in Figure 2a. As a result, we could conclude that the two systems above had similar laser transmission efficiency. The fluorescence intensity I 0 as a function of laser source power was measured in the confocal system (NA 2 = 0.9, objective) and the fiber system. As Figure 2b illustrated, the fluorescence intensity detected in the fiber system was about 1.87 times higher than that in the confocal system with the same laser source power. The sensitivity of the quantum magnetometer was inversely proportional to I 0 1/2 [25]. Therefore, the sensitivity performance of the fiber optical system was 27% higher than that of the confocal system.

Optical Interrogation of NV Centers in Diamond
The all-fiber NV center magnetometer was analyzed in terms of optical properties and sensitivity performance, compared with the traditional confocal system. Laser loss in the fiber system mainly occurred at AOM, fiber mating sleeves(FMS), and WMC. The diffraction efficiency of fiber-coupled AOM was 70%. Multi-mode fiber connectors provided low insertion loss, typically less than 0.5 dB, and the insertion loss from single-mode to multi-mode was negligible. The laser power transmitted to the diamond reached 64.2 mW when the laser source power was set to 100 mW, as measured by an optical power meter. In addition, the laser transmission efficiency of the traditional space light confocal system was about 62.3%, as shown in Figure 2a. As a result, we could conclude that the two systems above had similar laser transmission efficiency. The fluorescence intensity I0 as a function of laser source power was measured in the confocal system (NA2 = 0.9, objective) and the fiber system. As Figure 2b illustrated, the fluorescence intensity detected in the fiber system was about 1.87 times higher than that in the confocal system with the same laser source power. The sensitivity of the quantum magnetometer was inversely proportional to I0 1/2 [25]. Therefore, the sensitivity performance of the fiber optical system was 27% higher than that of the confocal system. The fluorescence intensity I0 as function of laser source power in confocal system and fiber-based system. The same micro-diamond was used to reduce the influence of diamond morphology. The fluorescence intensity is described using APD output voltage value.
As shown in Figure 3a, the diamond region (the red area) illuminated by the laser is a circular truncated cone. The diameter of the upper surface is a0 ≈ acore = 200 μm, and the divergence angle and acceptance angle are both θd = arcsin(NA1/nd) = 9.23°, where refractive index of diamond nd = 2.43 and the numerical aperture of Port2 NA1 = 0.39. The collecting efficiency of the optical fiber is estimated using geometric optics and the Fresnel reflection law [26]. The collecting efficiency of the fiber system is given by Equation (1). NV centers in micro-diamonds could be treated as a uniformly radiating point light source. The emitted fluorescence intensity distribution (Ip) of NV centers is 1/4π. The fluorescence intensity I 0 as function of laser source power in confocal system and fiber-based system. The same micro-diamond was used to reduce the influence of diamond morphology. The fluorescence intensity is described using APD output voltage value.
As shown in Figure 3a, the diamond region (the red area) illuminated by the laser is a circular truncated cone. The diameter of the upper surface is a 0 ≈ a core = 200 µm, and the divergence angle and acceptance angle are both θ d = arcsin(NA 1 /n d ) = 9.23 • , where refractive index of diamond n d = 2.43 and the numerical aperture of Port2 NA 1 = 0.39. The collecting efficiency of the optical fiber is estimated using geometric optics and the Fresnel reflection law [26]. The collecting efficiency of the fiber system is given by Equation (1). NV centers in micro-diamonds could be treated as a uniformly radiating point light source. The emitted fluorescence intensity distribution (I p ) of NV centers is 1/4π.
The collecting efficiency of the fiber system is estimated as 0.607% with a refractive index of the fiber core n core = 1.459. For the laser-excited region in the diamond, the laser power through the axial cross-section is consistent. In each tiny cylinder of the excited region with a thickness of dz, the power of the emitted fluorescence is (λ laser /λ fluo )δηN 0 P laser , where λ laser = 532 nm and the fluorescence wavelength range of NV center is 600-800 nm (Figure 3b). Here we used λ fluo = 700 nm to estimate the fluorescence intensity. δ ≈ 3 × 10 −17 cm 2 is the absorption cross section of type Ib diamond, η ≈ 0.1 is the PL quantum yield [27], N 0 = 1.08 × 10 18 /cm 3 is the NV center density, and P laser = 64.2 mW is the laser power exciting the diamond. The fluorescence power collected by the fiber is written as: where z d ≈ 550 µm is the micro-diamond dimension in the axial direction of the optical fiber. P f = 5 × 10 −2 mW is obtained by calculation with Equation (2). The optical loss from Port 2 of the WDM to the active area of the APD results in a transmission efficiency of about 87%. As a result, the fluorescence power detected by the APD is theoretically 4.35 × 10 −2 mW.
OR PEER REVIEW 5 of 12 The collecting efficiency of the fiber system is estimated as 0.607% with a refractive index of the fiber core ncore = 1.459. For the laser-excited region in the diamond, the laser power through the axial cross-section is consistent. In each tiny cylinder of the excited region with a thickness of dz, the power of the emitted fluorescence is (λlaser/λfluo)δηN0Plaser, where λlaser = 532 nm and the fluorescence wavelength range of NV center is 600-800 nm (Figure 3b). Here we used λfluo = 700 nm to estimate the fluorescence intensity. δ ≈ 3 × 10 −17 cm 2 is the absorption cross section of type Ib diamond, η ≈ 0.1 is the PL quantum yield [27], N0 = 1.08 × 10 18 /cm 3 is the NV center density, and Plaser = 64.2 mW is the laser power exciting the diamond. The fluorescence power collected by the fiber is written as: where zd ≈ 550 μm is the micro-diamond dimension in the axial direction of the optical fiber. Pf = 5 × 10 −2 mW is obtained by calculation with Equation (2). The optical loss from Port 2 of the WDM to the active area of the APD results in a transmission efficiency of about 87%. As a result, the fluorescence power detected by the APD is theoretically 4.35 × 10 −2 mW.
Experimentally, when the transimpedance gain was tuned to 2.15 × 10 4 v/w(700 nm), the APD detected fluorescence with an output voltage of 885.6 mV, which equates to 4.12 × 10 −2 mW. There is an acceptable deviation of 5.6% between the APD detected value and Experimentally, when the transimpedance gain was tuned to 2.15 × 10 4 v/w (700 nm), the APD detected fluorescence with an output voltage of 885.6 mV, which equates to 4.12 × 10 −2 mW. There is an acceptable deviation of 5.6% between the APD detected value and the theoretically estimated value. The following factors contribute to this difference. To begin with, there is a minor discrepancy between these parameters in Equation (2) and their actual values. Secondly, at the radial marginal region of the excited region, the fluorescence collecting efficiency is generally lower than 0.607%. In the confocal system, the excited area of diamond is the green part, as shown in Figure 4a. The collecting efficiency of the confocal system is given by Equation (3). Since the multi-mode fiber with a 200 µm fiber core serves as the pinhole, only the fluorescence radiating near the focus is collected and transmitted to the detector. The fluorescence acceptance angle is θ 0 = arcsin(NA 2 /n d ) = 21.74 • .
6 of 11 Nanomaterials 2023, 13, x FOR PEER REVIEW 6 of 12 The collecting efficiency of the confocal system is estimated as 2.937% with a refractive index of the fiber core nair = 1. An electronically controlled translation table controls the axial movement of a plane mirror to measure the axial optical sectioning curve of the confocal system, plotted in Figure 4b. The normalized axial optical sectioning curve indicates the axial optical transmission efficiency of fluorescence for every axial cross-section because defocus fluorescence signal is filtered by the pinhole [28]. The full width at halfmaximum(zc) is used to estimate the thickness at which diamond fluorescence is detected: zc = 62 μm. The fluorescence power collected by the objective is written as: According to Equations (2) and (4), we further deduced that: The objective lens of confocal system has a larger numerical aperture, which directly leads to a higher fluorescence collection efficiency. While the diamond volume integrated in the optical fiber system is much larger, and the final fluorescence intensity detected by the optical fiber system is higher. It can be concluded that there is a good consistency between the fluorescence collection theoretical calculation and experiment results for both the optical fiber system and the confocal system. This model could be used to estimate the optical interrogation of NV diamond. However, when the diamond size is close to or smaller than the fiber core diameter, parts of the laser irradiate out of the fiber and do not pump the diamond. The utilization rate of the laser decreases, and the model needs to be modified with the laser utilization rate.

Sensitivity Analysis
Pulsed ODMR experiments are implemented by measuring the normalized change of fluorescence intensity as a function of the microwave frequency Ω. When an external magnetic field B is applied to NV centers, Zeeman frequency splitting ∆Ω between ms = ±1 sublevels is induced. Since the fluorescence intensity from ms = ±1 state NV centers is weaker than ms = 0 state, when the microwave frequency Ω resonates with a transition between the ms = 0 and ms = ±1 sublevels, the fluorescence intensity decreases. Four pairs of resonance peaks Ωi± (i = 1, 2, 3, 4) would be generated corresponding to different projection (Bi) of the external magnetic field on four NV axes directions [29], as demonstrated in Equation (5). D and E are the zero-magnetic-field splitting and lattice stress parameter, The collecting efficiency of the confocal system is estimated as 2.937% with a refractive index of the fiber core n air = 1. An electronically controlled translation table controls the axial movement of a plane mirror to measure the axial optical sectioning curve of the confocal system, plotted in Figure 4b. The normalized axial optical sectioning curve indicates the axial optical transmission efficiency of fluorescence for every axial crosssection because defocus fluorescence signal is filtered by the pinhole [28]. The full width at half-maximum(z c ) is used to estimate the thickness at which diamond fluorescence is detected: z c = 62 µm. The fluorescence power collected by the objective is written as: According to Equations (2) and (4), we further deduced that: The objective lens of confocal system has a larger numerical aperture, which directly leads to a higher fluorescence collection efficiency. While the diamond volume integrated in the optical fiber system is much larger, and the final fluorescence intensity detected by the optical fiber system is higher. It can be concluded that there is a good consistency between the fluorescence collection theoretical calculation and experiment results for both the optical fiber system and the confocal system. This model could be used to estimate the optical interrogation of NV diamond. However, when the diamond size is close to or smaller than the fiber core diameter, parts of the laser irradiate out of the fiber and do not pump the diamond. The utilization rate of the laser decreases, and the model needs to be modified with the laser utilization rate.

Sensitivity Analysis
Pulsed ODMR experiments are implemented by measuring the normalized change of fluorescence intensity as a function of the microwave frequency Ω. When an external magnetic field B is applied to NV centers, Zeeman frequency splitting ∆Ω between m s = ±1 sublevels is induced. Since the fluorescence intensity from m s = ±1 state NV centers is weaker than m s = 0 state, when the microwave frequency Ω resonates with a transition between the m s = 0 and m s = ±1 sublevels, the fluorescence intensity decreases. Four pairs of resonance peaks Ω i± (i = 1, 2, 3, 4) would be generated corresponding to different projection (B i ) of the external magnetic field on four NV axes directions [29], as demonstrated in Equation (5). D and E are the zero-magnetic-field splitting and lattice stress parameter, γ = 2.8 MHz/Gauss is the gyromagnetic ratio. The magnitude of B can be measured directly through Ω i± according to Equation (6), in which ∆Ω i = Ω i+ − Ω i− .
The magnetometer sensitivity η is limited by photon shot-noise, which is estimated by Equation (7).
where ∆ν is the linewidth of the resonance peak in ODMR spectrum, C is the fluorescence intensity contrast, and R is the detected photon rate. The ODMR spectrum of the confocal system and fiber system were obtained with 20 Gauss DC magnetic field applied in the axial direction of the objective and fiber, which coincided with the <111> crystal direction of the diamond. As shown in Figure 5, the ODMR spectrum presented two pairs of resonance peaks, since the magnetic field along the <111> crystal direction has the same projection size on the other three NV axes. Lorentz fitting was performed on the resonance peaks in the box in Figure 5 to obtain the values of ∆ν 1 and C 1 for each system. For the confocal system, ∆ν 1 = 8.8 MHz and C 1 = 4.6% were obtained, while for the fiber system, ∆ν 2 = 10.2 MHz and C 2 = 3.2% were obtained. Wavelength of 700 nm was used to estimate the detected photon rate, R 1 = 0.78 × 10 14 Hz and R 2 = 1.45 × 10 14 Hz were obtained for confocal and fiber systems, respectively. Finally, the photon shot-noise limited sensitivity was 0.73 nT/Hz 1/2 for all-fiber magnetometer and 0.59 nT/Hz 1/2 for spatial optical confocal magnetometer.
The magnetometer sensitivity is limited by photon shot-noise, which is estimated by Equation (7).
where ∆ is the linewidth of the resonance peak in ODMR spectrum, C is the fluorescence intensity contrast, and R is the detected photon rate. The ODMR spectrum of the confocal system and fiber system were obtained with 20 Gauss DC magnetic field applied in the axial direction of the objective and fiber, which coincided with the <111> crystal direction of the diamond. As shown in Figure 5, the ODMR spectrum presented two pairs of resonance peaks, since the magnetic field along the <111> crystal direction has the same projection size on the other three NV axes. Lorentz fitting was performed on the resonance peaks in the box in Figure 5 to obtain the values of ∆ 1 and C1 for each system. For the confocal system, ∆ 1 = 8.8 MHz and C1 = 4.6% were obtained, while for the fiber system, ∆ 2 = 10.2 MHz and C2 = 3.2% were obtained. Wavelength of 700 nm was used to estimate the detected photon rate, R1 = 0.78 × 10 14 Hz and R2 = 1.45 × 10 14 Hz were obtained for confocal and fiber systems, respectively. Finally, the photon shot-noise limited sensitivity was 0.73 nT/Hz 1/2 for all-fiber magnetometer and 0.59 nT/Hz 1/2 for spatial optical confocal magnetometer.
Although, the all-fiber NV center magnetometer had higher fluorescence intensity. The all-fiber NV center magnetometer had a slightly lower sensitivity because of the wider linewidth and smaller contrast, compared with the confocal system. In the confocal system, laser was focused by the objective, resulting higher laser intensity distribution near the focal point where was the laser excited and fluorescence collected area of the diamond. The optical pumping rate of NV centers is proportional to the laser intensity. The ODMR spectra linewidth narrowing and contrast enhancement effect with increasing pumping rate has been discussed by Jensen et al. [30]. Compared with the fiber system, the linewidth of the confocal system deceased from 10.2 MHz to 8.8 MHz, and the contrast Although, the all-fiber NV center magnetometer had higher fluorescence intensity. The all-fiber NV center magnetometer had a slightly lower sensitivity because of the wider linewidth and smaller contrast, compared with the confocal system. In the confocal system, laser was focused by the objective, resulting higher laser intensity distribution near the focal point where was the laser excited and fluorescence collected area of the diamond. The optical pumping rate of NV centers is proportional to the laser intensity. The ODMR spectra linewidth narrowing and contrast enhancement effect with increasing pumping rate has been discussed by Jensen et al. [30]. Compared with the fiber system, the linewidth of the confocal system deceased from 10.2 MHz to 8.8 MHz, and the contrast boosted from 3.2% to 4.6%, which agreed with Jensen, K. et al. 's research result. Consequently, though a little sensitivity was sacrificed, the all-fiber magnetometer become robust and compact, could be applied for remote-magnetic measurement out of lab.

Vector Magnetic Field Measurement
There are challenges in defining the direction of the vector magnetic field in the lab frame. Firstly, the orientations of the four NV axes of the micro-diamond in lab frame are not easy to determine. Secondly, it is essential to distinguish these four NV axes and their corresponding resonance peaks and determine if B i is positive or negative [31].
Here, the 1 axis Helmholtz coil was introduced to solve the above-listed problems. The measurement steps of vector static magnetic field are illustrated in Figure 6. We set the angle between the auxiliary calibration magnetic field and the fiber to be 45 • (this angle could be set to be any value; here we set it to 45 • for the convenience of experiments and calculations). The C3v symmetry crystal lattice of diamond determines four possible orientations of NV center axes, as shown in Figure 7a. A Cartesian system of coordinates is built relative to fiber and the Helmholtz coil. One of the <111> crystal orientations in diamond, which is defined as axis 4, coincides with the axial direction of the optical fiber and is determined as z-axis. The x-, y-axes are introduced to set the auxiliary calibration magnetic field B a = (B 0 /2, B 0 /2, B 0 / √ 2) to facilitate testing and calibration. The projection of axes one, two, three and B a on the x-y plane is shown in Figure 7b. The axis one of NV centers is determined as its projection is the first axes that the x-axis encountered when it rotated counterclockwise in the x-y plane. We define the angle between axis one and the x-axis as θ (θ∈(0, 120 • )) in x-y plane. We define B NV = (B 1 , B 2 , B 3 , B 4 ) as the components in NV centers axes coordinate system, B lab = (B x , B y , B z ) as the components in lab-based Cartesian system. The transformation between the lab-based Cartesian system and the NV centers axes coordinate system is described as Equation (8). K is the transformation matrix formed by the unit vectors of the four NV centers axes.
boosted from 3.2% to 4.6%, which agreed with Jensen, K. et al. 's research result. Consequently, though a little sensitivity was sacrificed, the all-fiber magnetometer become robust and compact, could be applied for remote-magnetic measurement out of lab.

Vector Magnetic Field Measurement
There are challenges in defining the direction of the vector magnetic field in the lab frame. Firstly, the orientations of the four NV axes of the micro-diamond in lab frame are not easy to determine. Secondly, it is essential to distinguish these four NV axes and their corresponding resonance peaks and determine if Bi is positive or negative [31].
Here, the 1 axis Helmholtz coil was introduced to solve the above-listed problems. The measurement steps of vector static magnetic field are illustrated in Figure 6. We set the angle between the auxiliary calibration magnetic field and the fiber to be 45° (this angle could be set to be any value; here we set it to 45° for the convenience of experiments and calculations). The C3v symmetry crystal lattice of diamond determines four possible orientations of NV center axes, as shown in Figure 7a. A Cartesian system of coordinates is built relative to fiber and the Helmholtz coil. One of the <111> crystal orientations in diamond, which is defined as axis 4, coincides with the axial direction of the optical fiber and is determined as z-axis. The x-, y-axes are introduced to set the auxiliary calibration magnetic field Ba = (B0/2, B0/2, B0/√2) to facilitate testing and calibration. The projection of axes one, two, three and Ba on the x-y plane is shown in Figure 7b. The axis one of NV centers is determined as its projection is the first axes that the x-axis encountered when it rotated counterclockwise in the x-y plane. We define the angle between axis one and the x-axis as θ (θ∈(0, 120°)) in x-y plane. We define BNV = (B1, B2, B3, B4) as the components in NV centers axes coordinate system, Blab = (Bx, By, Bz) as the components in lab-based Cartesian system. The transformation between the lab-based Cartesian system and the NV centers axes coordinate system is described as Equation (8). K is the transformation matrix formed by the unit vectors of the four NV centers axes.   With an auxiliary calibration magnetic field B a generated by the Helmholtz coil with a 4A current, ODMR spectra were obtained, as Figure 8a shows. From the zero magnetic field measurements (shown as the inset in Figure 8a), D = 2871.5 MHz, E = 3.5 MHz were obtained. According to the result of the spectra fitting, we obtain differential resonance peaks With an auxiliary calibration magnetic field Ba generated by the Helmholtz coil with a 4A current, ODMR spectra were obtained, as Figure 8a shows. From the zero magnetic field measurements (shown as the inset in Figure 8a The ODMR spectrum of the static magnetic field to be measured Bm was shown in Figure 8b. To distinguish the NV axes and their corresponding resonance peaks, the auxiliary calibration magnetic field Ba was applied with Bm at the same time, and the ODMR spectrum of the magnetic field Ba + Bm was measured, as shown in Figure 8c. The resonance peaks in Figure 8c shifted compared to the spectrum in Figure 8b. The amount of shift for each peak is defined as dΩi = ∆Ωci − ∆Ωbi. We resolved the axes corresponding to each pair of resonance peaks through a numeric relationship |dΩi| = ∆Ωia. We can further determine if Bmi is positive or negative through ∆Ωia = ±dΩi. Finally, we could obtain BNVm = (1.   Figure 8. (a-c) are the ODMR spectra of calibration magnetic field Ba, measured magnetic field Bm, and superimposed magnetic field Ba + Bm, respectively. The bottom left inset of (a) shows the ODMR spectrum at zero magnetic field. The data in the figure are the frequency of resonance peaks Ωai, Ωbi, Ωci (i = 1, 2, 3, 4 are the four orientations of NV center axes).

Conclusions
We demonstrate an all-fiber magnetometer with a sensitivity of 0.73 nT/Hz 1/2 , which is portable and can be quickly built outside of a laboratory. Due to the good optical antiinterference of the fiber system, this sensor system could work in more complex working conditions which promotes the practical use of NV center sensors. This all-fiber magnetometer is a suitable candidate for endoscopic measurement of remote inaccessible places. This schematic could also be expanded to thermometer and magnetic field imaging. An optical interrogation model of NV centers in diamond was established, which was essen- Figure 8. (a-c) are the ODMR spectra of calibration magnetic field B a , measured magnetic field B m , and superimposed magnetic field B a + B m, respectively. The bottom left inset of (a) shows the ODMR spectrum at zero magnetic field. The data in the figure are the frequency of resonance peaks Ω ai , Ω bi , Ω ci (i = 1, 2, 3, 4 are the four orientations of NV center axes).
The ODMR spectrum of the static magnetic field to be measured B m was shown in Figure 8b. To distinguish the NV axes and their corresponding resonance peaks, the auxiliary calibration magnetic field B a was applied with B m at the same time, and the ODMR spectrum of the magnetic field B a + B m was measured, as shown in Figure 8c. The resonance peaks in Figure 8c shifted compared to the spectrum in Figure 8b. The amount of shift for each peak is defined as dΩ i = ∆Ω ci − ∆Ω bi . We resolved the axes corresponding to each pair of resonance peaks through a numeric relationship |dΩ i | = ∆Ω ia . We can further determine if B mi is positive or negative through ∆Ω ia = ±dΩ i . Finally, we could obtain B NVm = (1.19, 6.23, −15.61, 8.2) Gauss in the NV centers axis frame, corresponding to B labm = (−12.3, 6.7, 8.2) Gauss.

Conclusions
We demonstrate an all-fiber magnetometer with a sensitivity of 0.73 nT/Hz 1/2 , which is portable and can be quickly built outside of a laboratory. Due to the good optical anti-interference of the fiber system, this sensor system could work in more complex working conditions which promotes the practical use of NV center sensors. This all-fiber magnetometer is a suitable candidate for endoscopic measurement of remote inaccessible places. This schematic could also be expanded to thermometer and magnetic field imaging. An optical interrogation model of NV centers in diamond was established, which was essential to estimate the optical performance of the NV center system. A simplified rapid vector magnetic field measurement method is proposed for diamond particles with a one-axis Helmholtz coil, realizing µm-scale vector magnetic field detection.  Data Availability Statement: Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.